\(\int x^2 \cos ^3(a+b \log (c x^n)) \, dx\) [96]
Optimal result
Integrand size = 17, antiderivative size = 160 \[
\int x^2 \cos ^3\left (a+b \log \left (c x^n\right )\right ) \, dx=\frac {2 b^2 n^2 x^3 \cos \left (a+b \log \left (c x^n\right )\right )}{9+10 b^2 n^2+b^4 n^4}+\frac {x^3 \cos ^3\left (a+b \log \left (c x^n\right )\right )}{3 \left (1+b^2 n^2\right )}+\frac {2 b^3 n^3 x^3 \sin \left (a+b \log \left (c x^n\right )\right )}{3 \left (9+10 b^2 n^2+b^4 n^4\right )}+\frac {b n x^3 \cos ^2\left (a+b \log \left (c x^n\right )\right ) \sin \left (a+b \log \left (c x^n\right )\right )}{3 \left (1+b^2 n^2\right )}
\]
[Out]
2*b^2*n^2*x^3*cos(a+b*ln(c*x^n))/(b^4*n^4+10*b^2*n^2+9)+1/3*x^3*cos(a+b*ln(c*x^n))^3/(b^2*n^2+1)+2/3*b^3*n^3*x
^3*sin(a+b*ln(c*x^n))/(b^4*n^4+10*b^2*n^2+9)+1/3*b*n*x^3*cos(a+b*ln(c*x^n))^2*sin(a+b*ln(c*x^n))/(b^2*n^2+1)
Rubi [A] (verified)
Time = 0.07 (sec) , antiderivative size = 160, normalized size of antiderivative = 1.00, number
of steps used = 2, number of rules used = 2, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.118, Rules used = {4576, 4574}
\[
\int x^2 \cos ^3\left (a+b \log \left (c x^n\right )\right ) \, dx=\frac {x^3 \cos ^3\left (a+b \log \left (c x^n\right )\right )}{3 \left (b^2 n^2+1\right )}+\frac {b n x^3 \sin \left (a+b \log \left (c x^n\right )\right ) \cos ^2\left (a+b \log \left (c x^n\right )\right )}{3 \left (b^2 n^2+1\right )}+\frac {2 b^2 n^2 x^3 \cos \left (a+b \log \left (c x^n\right )\right )}{b^4 n^4+10 b^2 n^2+9}+\frac {2 b^3 n^3 x^3 \sin \left (a+b \log \left (c x^n\right )\right )}{3 \left (b^4 n^4+10 b^2 n^2+9\right )}
\]
[In]
Int[x^2*Cos[a + b*Log[c*x^n]]^3,x]
[Out]
(2*b^2*n^2*x^3*Cos[a + b*Log[c*x^n]])/(9 + 10*b^2*n^2 + b^4*n^4) + (x^3*Cos[a + b*Log[c*x^n]]^3)/(3*(1 + b^2*n
^2)) + (2*b^3*n^3*x^3*Sin[a + b*Log[c*x^n]])/(3*(9 + 10*b^2*n^2 + b^4*n^4)) + (b*n*x^3*Cos[a + b*Log[c*x^n]]^2
*Sin[a + b*Log[c*x^n]])/(3*(1 + b^2*n^2))
Rule 4574
Int[Cos[((a_.) + Log[(c_.)*(x_)^(n_.)]*(b_.))*(d_.)]*((e_.)*(x_))^(m_.), x_Symbol] :> Simp[(m + 1)*(e*x)^(m +
1)*(Cos[d*(a + b*Log[c*x^n])]/(b^2*d^2*e*n^2 + e*(m + 1)^2)), x] + Simp[b*d*n*(e*x)^(m + 1)*(Sin[d*(a + b*Log[
c*x^n])]/(b^2*d^2*e*n^2 + e*(m + 1)^2)), x] /; FreeQ[{a, b, c, d, e, m, n}, x] && NeQ[b^2*d^2*n^2 + (m + 1)^2,
0]
Rule 4576
Int[Cos[((a_.) + Log[(c_.)*(x_)^(n_.)]*(b_.))*(d_.)]^(p_)*((e_.)*(x_))^(m_.), x_Symbol] :> Simp[(m + 1)*(e*x)^
(m + 1)*(Cos[d*(a + b*Log[c*x^n])]^p/(b^2*d^2*e*n^2*p^2 + e*(m + 1)^2)), x] + (Dist[b^2*d^2*n^2*p*((p - 1)/(b^
2*d^2*n^2*p^2 + (m + 1)^2)), Int[(e*x)^m*Cos[d*(a + b*Log[c*x^n])]^(p - 2), x], x] + Simp[b*d*n*p*(e*x)^(m + 1
)*Sin[d*(a + b*Log[c*x^n])]*(Cos[d*(a + b*Log[c*x^n])]^(p - 1)/(b^2*d^2*e*n^2*p^2 + e*(m + 1)^2)), x]) /; Free
Q[{a, b, c, d, e, m, n}, x] && IGtQ[p, 1] && NeQ[b^2*d^2*n^2*p^2 + (m + 1)^2, 0]
Rubi steps \begin{align*}
\text {integral}& = \frac {x^3 \cos ^3\left (a+b \log \left (c x^n\right )\right )}{3 \left (1+b^2 n^2\right )}+\frac {b n x^3 \cos ^2\left (a+b \log \left (c x^n\right )\right ) \sin \left (a+b \log \left (c x^n\right )\right )}{3 \left (1+b^2 n^2\right )}+\frac {\left (2 b^2 n^2\right ) \int x^2 \cos \left (a+b \log \left (c x^n\right )\right ) \, dx}{3 \left (1+b^2 n^2\right )} \\ & = \frac {2 b^2 n^2 x^3 \cos \left (a+b \log \left (c x^n\right )\right )}{9+10 b^2 n^2+b^4 n^4}+\frac {x^3 \cos ^3\left (a+b \log \left (c x^n\right )\right )}{3 \left (1+b^2 n^2\right )}+\frac {2 b^3 n^3 x^3 \sin \left (a+b \log \left (c x^n\right )\right )}{3 \left (9+10 b^2 n^2+b^4 n^4\right )}+\frac {b n x^3 \cos ^2\left (a+b \log \left (c x^n\right )\right ) \sin \left (a+b \log \left (c x^n\right )\right )}{3 \left (1+b^2 n^2\right )} \\
\end{align*}
Mathematica [A] (verified)
Time = 0.38 (sec) , antiderivative size = 120, normalized size of antiderivative = 0.75
\[
\int x^2 \cos ^3\left (a+b \log \left (c x^n\right )\right ) \, dx=\frac {x^3 \left (27 \left (1+b^2 n^2\right ) \cos \left (a+b \log \left (c x^n\right )\right )+\left (9+b^2 n^2\right ) \cos \left (3 \left (a+b \log \left (c x^n\right )\right )\right )+2 b n \left (9+5 b^2 n^2+\left (9+b^2 n^2\right ) \cos \left (2 \left (a+b \log \left (c x^n\right )\right )\right )\right ) \sin \left (a+b \log \left (c x^n\right )\right )\right )}{12 \left (9+10 b^2 n^2+b^4 n^4\right )}
\]
[In]
Integrate[x^2*Cos[a + b*Log[c*x^n]]^3,x]
[Out]
(x^3*(27*(1 + b^2*n^2)*Cos[a + b*Log[c*x^n]] + (9 + b^2*n^2)*Cos[3*(a + b*Log[c*x^n])] + 2*b*n*(9 + 5*b^2*n^2
+ (9 + b^2*n^2)*Cos[2*(a + b*Log[c*x^n])])*Sin[a + b*Log[c*x^n]]))/(12*(9 + 10*b^2*n^2 + b^4*n^4))
Maple [A] (verified)
Time = 7.31 (sec) , antiderivative size = 117, normalized size of antiderivative = 0.73
| | |
| method | result | size |
| | |
| parallelrisch |
\(\frac {9 x^{3} \left (\left (\frac {b^{2} n^{2}}{9}+1\right ) \cos \left (3 b \ln \left (c \,x^{n}\right )+3 a \right )+\left (\frac {1}{9} b^{3} n^{3}+b n \right ) \sin \left (3 b \ln \left (c \,x^{n}\right )+3 a \right )+\left (b^{2} n^{2}+1\right ) \left (\sin \left (a +b \ln \left (c \,x^{n}\right )\right ) b n +3 \cos \left (a +b \ln \left (c \,x^{n}\right )\right )\right )\right )}{12 b^{4} n^{4}+120 b^{2} n^{2}+108}\) |
\(117\) |
| | |
|
|
|
|
[In]
int(x^2*cos(a+b*ln(c*x^n))^3,x,method=_RETURNVERBOSE)
[Out]
9*x^3*((1/9*b^2*n^2+1)*cos(3*b*ln(c*x^n)+3*a)+(1/9*b^3*n^3+b*n)*sin(3*b*ln(c*x^n)+3*a)+(b^2*n^2+1)*(sin(a+b*ln
(c*x^n))*b*n+3*cos(a+b*ln(c*x^n))))/(12*b^4*n^4+120*b^2*n^2+108)
Fricas [A] (verification not implemented)
none
Time = 0.25 (sec) , antiderivative size = 127, normalized size of antiderivative = 0.79
\[
\int x^2 \cos ^3\left (a+b \log \left (c x^n\right )\right ) \, dx=\frac {6 \, b^{2} n^{2} x^{3} \cos \left (b n \log \left (x\right ) + b \log \left (c\right ) + a\right ) + {\left (b^{2} n^{2} + 9\right )} x^{3} \cos \left (b n \log \left (x\right ) + b \log \left (c\right ) + a\right )^{3} + {\left (2 \, b^{3} n^{3} x^{3} + {\left (b^{3} n^{3} + 9 \, b n\right )} x^{3} \cos \left (b n \log \left (x\right ) + b \log \left (c\right ) + a\right )^{2}\right )} \sin \left (b n \log \left (x\right ) + b \log \left (c\right ) + a\right )}{3 \, {\left (b^{4} n^{4} + 10 \, b^{2} n^{2} + 9\right )}}
\]
[In]
integrate(x^2*cos(a+b*log(c*x^n))^3,x, algorithm="fricas")
[Out]
1/3*(6*b^2*n^2*x^3*cos(b*n*log(x) + b*log(c) + a) + (b^2*n^2 + 9)*x^3*cos(b*n*log(x) + b*log(c) + a)^3 + (2*b^
3*n^3*x^3 + (b^3*n^3 + 9*b*n)*x^3*cos(b*n*log(x) + b*log(c) + a)^2)*sin(b*n*log(x) + b*log(c) + a))/(b^4*n^4 +
10*b^2*n^2 + 9)
Sympy [F]
\[
\int x^2 \cos ^3\left (a+b \log \left (c x^n\right )\right ) \, dx=\begin {cases} \int x^{2} \cos ^{3}{\left (a - \frac {3 i \log {\left (c x^{n} \right )}}{n} \right )}\, dx & \text {for}\: b = - \frac {3 i}{n} \\\int x^{2} \cos ^{3}{\left (a - \frac {i \log {\left (c x^{n} \right )}}{n} \right )}\, dx & \text {for}\: b = - \frac {i}{n} \\\int x^{2} \cos ^{3}{\left (a + \frac {i \log {\left (c x^{n} \right )}}{n} \right )}\, dx & \text {for}\: b = \frac {i}{n} \\\int x^{2} \cos ^{3}{\left (a + \frac {3 i \log {\left (c x^{n} \right )}}{n} \right )}\, dx & \text {for}\: b = \frac {3 i}{n} \\\frac {2 b^{3} n^{3} x^{3} \sin ^{3}{\left (a + b \log {\left (c x^{n} \right )} \right )}}{3 b^{4} n^{4} + 30 b^{2} n^{2} + 27} + \frac {3 b^{3} n^{3} x^{3} \sin {\left (a + b \log {\left (c x^{n} \right )} \right )} \cos ^{2}{\left (a + b \log {\left (c x^{n} \right )} \right )}}{3 b^{4} n^{4} + 30 b^{2} n^{2} + 27} + \frac {6 b^{2} n^{2} x^{3} \sin ^{2}{\left (a + b \log {\left (c x^{n} \right )} \right )} \cos {\left (a + b \log {\left (c x^{n} \right )} \right )}}{3 b^{4} n^{4} + 30 b^{2} n^{2} + 27} + \frac {7 b^{2} n^{2} x^{3} \cos ^{3}{\left (a + b \log {\left (c x^{n} \right )} \right )}}{3 b^{4} n^{4} + 30 b^{2} n^{2} + 27} + \frac {9 b n x^{3} \sin {\left (a + b \log {\left (c x^{n} \right )} \right )} \cos ^{2}{\left (a + b \log {\left (c x^{n} \right )} \right )}}{3 b^{4} n^{4} + 30 b^{2} n^{2} + 27} + \frac {9 x^{3} \cos ^{3}{\left (a + b \log {\left (c x^{n} \right )} \right )}}{3 b^{4} n^{4} + 30 b^{2} n^{2} + 27} & \text {otherwise} \end {cases}
\]
[In]
integrate(x**2*cos(a+b*ln(c*x**n))**3,x)
[Out]
Piecewise((Integral(x**2*cos(a - 3*I*log(c*x**n)/n)**3, x), Eq(b, -3*I/n)), (Integral(x**2*cos(a - I*log(c*x**
n)/n)**3, x), Eq(b, -I/n)), (Integral(x**2*cos(a + I*log(c*x**n)/n)**3, x), Eq(b, I/n)), (Integral(x**2*cos(a
+ 3*I*log(c*x**n)/n)**3, x), Eq(b, 3*I/n)), (2*b**3*n**3*x**3*sin(a + b*log(c*x**n))**3/(3*b**4*n**4 + 30*b**2
*n**2 + 27) + 3*b**3*n**3*x**3*sin(a + b*log(c*x**n))*cos(a + b*log(c*x**n))**2/(3*b**4*n**4 + 30*b**2*n**2 +
27) + 6*b**2*n**2*x**3*sin(a + b*log(c*x**n))**2*cos(a + b*log(c*x**n))/(3*b**4*n**4 + 30*b**2*n**2 + 27) + 7*
b**2*n**2*x**3*cos(a + b*log(c*x**n))**3/(3*b**4*n**4 + 30*b**2*n**2 + 27) + 9*b*n*x**3*sin(a + b*log(c*x**n))
*cos(a + b*log(c*x**n))**2/(3*b**4*n**4 + 30*b**2*n**2 + 27) + 9*x**3*cos(a + b*log(c*x**n))**3/(3*b**4*n**4 +
30*b**2*n**2 + 27), True))
Maxima [B] (verification not implemented)
Leaf count of result is larger than twice the leaf count of optimal. 1007 vs. \(2 (154) = 308\).
Time = 0.26 (sec) , antiderivative size = 1007, normalized size of antiderivative = 6.29
\[
\int x^2 \cos ^3\left (a+b \log \left (c x^n\right )\right ) \, dx=\text {Too large to display}
\]
[In]
integrate(x^2*cos(a+b*log(c*x^n))^3,x, algorithm="maxima")
[Out]
1/24*(((b^3*cos(3*b*log(c))*sin(6*b*log(c)) - b^3*cos(6*b*log(c))*sin(3*b*log(c)) + b^3*sin(3*b*log(c)))*n^3 +
(b^2*cos(6*b*log(c))*cos(3*b*log(c)) + b^2*sin(6*b*log(c))*sin(3*b*log(c)) + b^2*cos(3*b*log(c)))*n^2 + 9*(b*
cos(3*b*log(c))*sin(6*b*log(c)) - b*cos(6*b*log(c))*sin(3*b*log(c)) + b*sin(3*b*log(c)))*n + 9*cos(6*b*log(c))
*cos(3*b*log(c)) + 9*sin(6*b*log(c))*sin(3*b*log(c)) + 9*cos(3*b*log(c)))*x^3*cos(3*b*log(x^n) + 3*a) + 9*((b^
3*cos(3*b*log(c))*sin(4*b*log(c)) - b^3*cos(4*b*log(c))*sin(3*b*log(c)) + b^3*cos(2*b*log(c))*sin(3*b*log(c))
- b^3*cos(3*b*log(c))*sin(2*b*log(c)))*n^3 + 3*(b^2*cos(4*b*log(c))*cos(3*b*log(c)) + b^2*cos(3*b*log(c))*cos(
2*b*log(c)) + b^2*sin(4*b*log(c))*sin(3*b*log(c)) + b^2*sin(3*b*log(c))*sin(2*b*log(c)))*n^2 + (b*cos(3*b*log(
c))*sin(4*b*log(c)) - b*cos(4*b*log(c))*sin(3*b*log(c)) + b*cos(2*b*log(c))*sin(3*b*log(c)) - b*cos(3*b*log(c)
)*sin(2*b*log(c)))*n + 3*cos(4*b*log(c))*cos(3*b*log(c)) + 3*cos(3*b*log(c))*cos(2*b*log(c)) + 3*sin(4*b*log(c
))*sin(3*b*log(c)) + 3*sin(3*b*log(c))*sin(2*b*log(c)))*x^3*cos(b*log(x^n) + a) + ((b^3*cos(6*b*log(c))*cos(3*
b*log(c)) + b^3*sin(6*b*log(c))*sin(3*b*log(c)) + b^3*cos(3*b*log(c)))*n^3 - (b^2*cos(3*b*log(c))*sin(6*b*log(
c)) - b^2*cos(6*b*log(c))*sin(3*b*log(c)) + b^2*sin(3*b*log(c)))*n^2 + 9*(b*cos(6*b*log(c))*cos(3*b*log(c)) +
b*sin(6*b*log(c))*sin(3*b*log(c)) + b*cos(3*b*log(c)))*n - 9*cos(3*b*log(c))*sin(6*b*log(c)) + 9*cos(6*b*log(c
))*sin(3*b*log(c)) - 9*sin(3*b*log(c)))*x^3*sin(3*b*log(x^n) + 3*a) + 9*((b^3*cos(4*b*log(c))*cos(3*b*log(c))
+ b^3*cos(3*b*log(c))*cos(2*b*log(c)) + b^3*sin(4*b*log(c))*sin(3*b*log(c)) + b^3*sin(3*b*log(c))*sin(2*b*log(
c)))*n^3 - 3*(b^2*cos(3*b*log(c))*sin(4*b*log(c)) - b^2*cos(4*b*log(c))*sin(3*b*log(c)) + b^2*cos(2*b*log(c))*
sin(3*b*log(c)) - b^2*cos(3*b*log(c))*sin(2*b*log(c)))*n^2 + (b*cos(4*b*log(c))*cos(3*b*log(c)) + b*cos(3*b*lo
g(c))*cos(2*b*log(c)) + b*sin(4*b*log(c))*sin(3*b*log(c)) + b*sin(3*b*log(c))*sin(2*b*log(c)))*n - 3*cos(3*b*l
og(c))*sin(4*b*log(c)) + 3*cos(4*b*log(c))*sin(3*b*log(c)) - 3*cos(2*b*log(c))*sin(3*b*log(c)) + 3*cos(3*b*log
(c))*sin(2*b*log(c)))*x^3*sin(b*log(x^n) + a))/((b^4*cos(3*b*log(c))^2 + b^4*sin(3*b*log(c))^2)*n^4 + 10*(b^2*
cos(3*b*log(c))^2 + b^2*sin(3*b*log(c))^2)*n^2 + 9*cos(3*b*log(c))^2 + 9*sin(3*b*log(c))^2)
Giac [B] (verification not implemented)
Leaf count of result is larger than twice the leaf count of optimal. 18053 vs. \(2 (154) = 308\).
Time = 1.46 (sec) , antiderivative size = 18053, normalized size of antiderivative = 112.83
\[
\int x^2 \cos ^3\left (a+b \log \left (c x^n\right )\right ) \, dx=\text {Too large to display}
\]
[In]
integrate(x^2*cos(a+b*log(c*x^n))^3,x, algorithm="giac")
[Out]
-1/24*(18*b^3*n^3*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x))
+ 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a) + 18*b^3*n^3*x^3
*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))
^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a) + 2*b^3*n^3*x^3*e^(3/2*pi*b*n*sgn(x)
- 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs
(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)*tan(1/2*a)^2 + 2*b^3*n^3*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*p
i*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(
c)))^2*tan(3/2*a)*tan(1/2*a)^2 + 18*b^3*n^3*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b
)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(3/2*a)^2*tan
(1/2*a)^2 + 18*b^3*n^3*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*b*n*log(ab
s(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(3/2*a)^2*tan(1/2*a)^2 + 2*b^3*n^
3*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(
c)))*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2 + 2*b^3*n^3*x^3*e^(-3/2*pi*b*n*s
gn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(1/2*b*n*log(
abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2 - b^2*n^2*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2
*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(ab
s(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2 - 27*b^2*n^2*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*
pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)
^2*tan(1/2*a)^2 - 27*b^2*n^2*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*b*n*
log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2 -
b^2*n^2*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*l
og(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2 + 2*b^3*n^3*x^3*e^(3/2*
pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/
2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a) + 2*b^3*n^3*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*p
i*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(
c)))^2*tan(3/2*a) - 18*b^3*n^3*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*b*n
*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(3/2*a)^2 - 18*b^3*n^3*x^3
*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))
^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(3/2*a)^2 + 2*b^3*n^3*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n
+ 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(1/2*b*n*log(abs(x)) + 1/2*b*lo
g(abs(c)))^2*tan(3/2*a)^2 + 2*b^3*n^3*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan
(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2 + 18*b^3
*n^3*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(a
bs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(1/2*a) + 18*b^3*n^3*x^3*e^(-1/2*pi*b*n*sgn(x) + 1
/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x))
+ 1/2*b*log(abs(c)))^2*tan(1/2*a) - 18*b^3*n^3*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*
pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a) - 18*b^3*n^3*x^3*e^(-1/2*pi*b*n*s
gn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(3/2*a)^2*t
an(1/2*a) + 18*b^3*n^3*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(1/2*b*n*log(abs
(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a) + 18*b^3*n^3*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*
pi*b*sgn(c) + 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a) + 18*b^3*n^3*x^
3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))
^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(1/2*a)^2 + 18*b^3*n^3*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b
*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*
b*log(abs(c)))*tan(1/2*a)^2 - 2*b^3*n^3*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*ta
n(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(1/2*a)^2 - 2*b^3
*n^3*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(
abs(c)))*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(1/2*a)^2 + 2*b^3*n^3*x^3*e^(3/2*pi*b*n*sgn(x) - 3/
2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(3/2*a)*tan(1/2*a)^2
+ 2*b^3*n^3*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2
*b*log(abs(c)))^2*tan(3/2*a)*tan(1/2*a)^2 - 2*b^3*n^3*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c)
- 3/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)*tan(1/2*a)^2 - 2*b^3*n^3*x^3*e^(-3/2*pi*
b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a
)*tan(1/2*a)^2 + 2*b^3*n^3*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log
(abs(x)) + 3/2*b*log(abs(c)))*tan(3/2*a)^2*tan(1/2*a)^2 + 2*b^3*n^3*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3
/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(3/2*a)^2*tan(1/2*a)^2 + 18*b^3*n^3
*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c
)))*tan(3/2*a)^2*tan(1/2*a)^2 + 18*b^3*n^3*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b
)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(3/2*a)^2*tan(1/2*a)^2 - b^2*n^2*x^3*e^(3/2*pi*b*n*sgn(x) -
3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)
) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2 + 27*b^2*n^2*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1
/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2
*a)^2 + 27*b^2*n^2*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*b*n*log(abs(x)
) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2 - b^2*n^2*x^3*e^(-3/2*pi*
b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b
*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2 + 108*b^2*n^2*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*p
i*b*sgn(c) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(
c)))*tan(3/2*a)^2*tan(1/2*a) + 108*b^2*n^2*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b
)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(3/2*a)^2*tan
(1/2*a) + b^2*n^2*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x))
+ 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(1/2*a)^2 - 27*b^2*n^2*x^3*e^(1/2*pi*
b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b
*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(1/2*a)^2 - 27*b^2*n^2*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*p
i*b*sgn(c) + 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(
c)))^2*tan(1/2*a)^2 + b^2*n^2*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n
*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(1/2*a)^2 + 4*b^2*n^2*x^
3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))
*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)*tan(1/2*a)^2 + 4*b^2*n^2*x^3*e^(-3/2*pi*b*n*sgn(x)
+ 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(1/2*b*n*log(abs(x)
) + 1/2*b*log(abs(c)))^2*tan(3/2*a)*tan(1/2*a)^2 - b^2*n^2*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sg
n(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2 + 27*b^2*n^2*x^3*e^(
1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*ta
n(3/2*a)^2*tan(1/2*a)^2 + 27*b^2*n^2*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(
3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2 - b^2*n^2*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2
*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2
+ b^2*n^2*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b
*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2 - 27*b^2*n^2*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c)
- 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2 - 27*b^2*n^2*x^3*e^(-1/2
*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3
/2*a)^2*tan(1/2*a)^2 + b^2*n^2*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(1/2*b*
n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2 - 18*b^3*n^3*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*
b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2
*b*log(abs(c))) - 18*b^3*n^3*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*b*n*
log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c))) - 2*b^3*n^3*x^3*e^(3/2*pi*b*n*
sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(1/2*b*n*log
(abs(x)) + 1/2*b*log(abs(c)))^2 - 2*b^3*n^3*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*
b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2 + 2*b^3*n^3*x^3
*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^
2*tan(3/2*a) + 2*b^3*n^3*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(
abs(x)) + 3/2*b*log(abs(c)))^2*tan(3/2*a) - 2*b^3*n^3*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c)
- 3/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a) - 2*b^3*n^3*x^3*e^(-3/2*pi*b*n*sgn(x) +
3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a) + 2*b^3*n^3
*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c
)))*tan(3/2*a)^2 + 2*b^3*n^3*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*
log(abs(x)) + 3/2*b*log(abs(c)))*tan(3/2*a)^2 - 18*b^3*n^3*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sg
n(c) - 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(3/2*a)^2 - 18*b^3*n^3*x^3*e^(-1/2*pi*b*n*sgn
(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(3/2*a)^2 - 18*
b^3*n^3*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*lo
g(abs(c)))^2*tan(1/2*a) - 18*b^3*n^3*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(
3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*a) + 18*b^3*n^3*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2
*pi*b*sgn(c) - 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(1/2*a) + 18*b^3*n^3*x^3*e^(-1/2*pi
*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(1/2*
a) - 18*b^3*n^3*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*a)^2*tan(1/2*a) -
18*b^3*n^3*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*a)^2*tan(1/2*a) + 18*b
*n*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs
(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a) + 18*b*n*x^3*e^(-1/2*pi*b*n*sgn
(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(
abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a) - 2*b^3*n^3*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2
*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(1/2*a)^2 - 2*b^3*n^3*x^3*e^(-3/2*pi*
b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(1/2*a)^
2 + 18*b^3*n^3*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1
/2*b*log(abs(c)))*tan(1/2*a)^2 + 18*b^3*n^3*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*
b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(1/2*a)^2 - 2*b^3*n^3*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n
+ 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*a)*tan(1/2*a)^2 - 2*b^3*n^3*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/
2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*a)*tan(1/2*a)^2 + 18*b*n*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sg
n(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2
*tan(3/2*a)*tan(1/2*a)^2 + 18*b*n*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2
*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)*tan(1/2*a)^2
+ 18*b*n*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*
log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(3/2*a)^2*tan(1/2*a)^2 + 18*b*n*x^3*e^(-1/2*pi*
b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b
*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(3/2*a)^2*tan(1/2*a)^2 + 18*b*n*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n +
3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(
abs(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2 + 18*b*n*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*p
i*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2*
tan(1/2*a)^2 + b^2*n^2*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs
(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2 + 27*b^2*n^2*x^3*e^(1/2*pi*b*n*sgn(
x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(a
bs(x)) + 1/2*b*log(abs(c)))^2 + 27*b^2*n^2*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b
)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2 + b^2*n^2*x^3*
e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^
2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2 + 4*b^2*n^2*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*
sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2
*tan(3/2*a) + 4*b^2*n^2*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(a
bs(x)) + 3/2*b*log(abs(c)))*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a) - b^2*n^2*x^3*e^(3/2*pi*
b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(3/2*a
)^2 - 27*b^2*n^2*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) +
3/2*b*log(abs(c)))^2*tan(3/2*a)^2 - 27*b^2*n^2*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2
*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(3/2*a)^2 - b^2*n^2*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*p
i*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(3/2*a)^2 + b^2*n^2*x^3*
e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2
*tan(3/2*a)^2 + 27*b^2*n^2*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(1/2*b*n*log
(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2 + 27*b^2*n^2*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sg
n(c) + 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2 + b^2*n^2*x^3*e^(-3/2*pi*b*n*sgn(
x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2 + 10
8*b^2*n^2*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*
log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(1/2*a) + 108*b^2*n^2*x^3*e^(-1/2*pi*b*n*sgn(x)
+ 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs
(x)) + 1/2*b*log(abs(c)))*tan(1/2*a) + 108*b^2*n^2*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1
/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(3/2*a)^2*tan(1/2*a) + 108*b^2*n^2*x^3*e^(-1/2*pi*b*n
*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(3/2*a)^2*t
an(1/2*a) + b^2*n^2*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)
) + 3/2*b*log(abs(c)))^2*tan(1/2*a)^2 + 27*b^2*n^2*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1
/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*a)^2 + 27*b^2*n^2*x^3*e^(-1/2*pi*b*n*sgn(x) +
1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*a)^2 + b^2*n^2
*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(
c)))^2*tan(1/2*a)^2 - b^2*n^2*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(1/2*b*n*
log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(1/2*a)^2 - 27*b^2*n^2*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*
sgn(c) - 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(1/2*a)^2 - 27*b^2*n^2*x^3*e^(-1/2*pi*b*n
*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(1/2*a)^2
- b^2*n^2*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*
b*log(abs(c)))^2*tan(1/2*a)^2 + 4*b^2*n^2*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*
tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(3/2*a)*tan(1/2*a)^2 + 4*b^2*n^2*x^3*e^(-3/2*pi*b*n*sgn(x) + 3
/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(3/2*a)*tan(1/2*a)^2 +
b^2*n^2*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*a)^2*tan(1/2*a)^2 + 27*b^
2*n^2*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*a)^2*tan(1/2*a)^2 + 27*b^2*n
^2*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*a)^2*tan(1/2*a)^2 + b^2*n^2*x^
3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*a)^2*tan(1/2*a)^2 - 9*x^3*e^(3/2*pi
*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*
b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2 - 27*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n +
1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log
(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2 - 27*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b
)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2*t
an(1/2*a)^2 - 9*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) +
3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2 - 2*b^3*n^3*x^3
*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))
- 2*b^3*n^3*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2
*b*log(abs(c))) - 18*b^3*n^3*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(1/2*b*n*l
og(abs(x)) + 1/2*b*log(abs(c))) - 18*b^3*n^3*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi
*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c))) - 2*b^3*n^3*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*
sgn(c) - 3/2*pi*b)*tan(3/2*a) - 2*b^3*n^3*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)
*tan(3/2*a) + 18*b*n*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x
)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a) + 18*b*n*x^3*e^(-3/2*pi*b*
n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n
*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a) - 18*b*n*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c
) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(
3/2*a)^2 - 18*b*n*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*b*n*log(abs(x))
+ 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(3/2*a)^2 + 18*b*n*x^3*e^(3/2*pi*b*n*s
gn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(1/2*b*n*log(
abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2 + 18*b*n*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c)
+ 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/
2*a)^2 - 18*b^3*n^3*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(1/2*a) - 18*b^3*n^
3*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(1/2*a) + 18*b*n*x^3*e^(1/2*pi*b*n*s
gn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*lo
g(abs(x)) + 1/2*b*log(abs(c)))^2*tan(1/2*a) + 18*b*n*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c)
+ 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(
1/2*a) - 18*b*n*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) +
3/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a) - 18*b*n*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c)
+ 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a) + 18*b*n*x^3*e^(1/2*pi*b*n
*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2
*tan(1/2*a) + 18*b*n*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(1/2*b*n*log(abs(
x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a) + 18*b*n*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*s
gn(c) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*
tan(1/2*a)^2 + 18*b*n*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*b*n*log(abs
(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(1/2*a)^2 - 18*b*n*x^3*e^(3/2*pi*b
*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(1/2*b*n*
log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(1/2*a)^2 - 18*b*n*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn
(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*ta
n(1/2*a)^2 + 18*b*n*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)
) + 3/2*b*log(abs(c)))^2*tan(3/2*a)*tan(1/2*a)^2 + 18*b*n*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sg
n(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(3/2*a)*tan(1/2*a)^2 - 18*b*n*x^3*e^(3/2*pi
*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*
a)*tan(1/2*a)^2 - 18*b*n*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(1/2*b*n*log(
abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)*tan(1/2*a)^2 + 18*b*n*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi
*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(3/2*a)^2*tan(1/2*a)^2 + 18*b*n*x^3*e^(-
3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(
3/2*a)^2*tan(1/2*a)^2 + 18*b*n*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(1/2*b*n
*log(abs(x)) + 1/2*b*log(abs(c)))*tan(3/2*a)^2*tan(1/2*a)^2 + 18*b*n*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n -
1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(3/2*a)^2*tan(1/2*a)^2 + b^2*n^2*x
^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c))
)^2 - 27*b^2*n^2*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) +
3/2*b*log(abs(c)))^2 - 27*b^2*n^2*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/
2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2 + b^2*n^2*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) +
3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2 - b^2*n^2*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2
*pi*b*sgn(c) - 3/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2 + 27*b^2*n^2*x^3*e^(1/2*pi*b*n*sgn(x)
- 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2 + 27*b^2*n^2*x^3*e^(
-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2 -
b^2*n^2*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*
log(abs(c)))^2 + 4*b^2*n^2*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log
(abs(x)) + 3/2*b*log(abs(c)))*tan(3/2*a) + 4*b^2*n^2*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c)
+ 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(3/2*a) + b^2*n^2*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*p
i*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*a)^2 - 27*b^2*n^2*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b
*sgn(c) - 1/2*pi*b)*tan(3/2*a)^2 - 27*b^2*n^2*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*p
i*b)*tan(3/2*a)^2 + b^2*n^2*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*a)^2
- 9*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(ab
s(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2 + 27*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*
n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b
*log(abs(c)))^2*tan(3/2*a)^2 + 27*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2
*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2 - 9*x^3*e^
(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*
tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2 + 108*b^2*n^2*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*
n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(1/2*a) + 108*b^2*n^2*x^3*e^(-
1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(
1/2*a) + 108*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2
*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(3/2*a)^2*tan(1/2*a) + 108*x^3*e^(-1/2*pi*b*
n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n
*log(abs(x)) + 1/2*b*log(abs(c)))*tan(3/2*a)^2*tan(1/2*a) - b^2*n^2*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/
2*pi*b*sgn(c) - 3/2*pi*b)*tan(1/2*a)^2 + 27*b^2*n^2*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) -
1/2*pi*b)*tan(1/2*a)^2 + 27*b^2*n^2*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(1
/2*a)^2 - b^2*n^2*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(1/2*a)^2 + 9*x^3*e^
(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*t
an(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(1/2*a)^2 - 27*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi
*b*sgn(c) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c
)))^2*tan(1/2*a)^2 - 27*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*b*n*log(a
bs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(1/2*a)^2 + 9*x^3*e^(-3/2*pi*b
*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*
n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(1/2*a)^2 + 36*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c)
- 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3
/2*a)*tan(1/2*a)^2 + 36*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(a
bs(x)) + 3/2*b*log(abs(c)))*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)*tan(1/2*a)^2 - 9*x^3*e^(
3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*ta
n(3/2*a)^2*tan(1/2*a)^2 + 27*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*b*n*l
og(abs(x)) + 3/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2 + 27*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*
pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2 - 9*x^3*e^(-3
/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan
(3/2*a)^2*tan(1/2*a)^2 + 9*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(1/2*b*n*log
(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2 - 27*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*
b*sgn(c) - 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2 - 27*x^3*e^(-1/2
*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3
/2*a)^2*tan(1/2*a)^2 + 9*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(1/2*b*n*log(
abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2 - 18*b*n*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*
pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs
(c))) - 18*b*n*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) +
3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c))) - 18*b*n*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b
*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(1/2*b*n*log(abs(x)) + 1/2*b*
log(abs(c)))^2 - 18*b*n*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(a
bs(x)) + 3/2*b*log(abs(c)))*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2 + 18*b*n*x^3*e^(3/2*pi*b*n*sgn(x) -
3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(3/2*a) + 18*b*n*x
^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)
))^2*tan(3/2*a) - 18*b*n*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(1/2*b*n*log(a
bs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a) - 18*b*n*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3
/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a) + 18*b*n*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*
b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(3/2*a)^2 + 18*b*n*x^3*e^(-3
/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(3
/2*a)^2 - 18*b*n*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) +
1/2*b*log(abs(c)))*tan(3/2*a)^2 - 18*b*n*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)
*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(3/2*a)^2 - 18*b*n*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/
2*pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*a) - 18*b*n*x^3*e^(-1/2*pi*b*
n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*a)
+ 18*b*n*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*l
og(abs(c)))^2*tan(1/2*a) + 18*b*n*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(1/2
*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(1/2*a) - 18*b*n*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*s
gn(c) - 1/2*pi*b)*tan(3/2*a)^2*tan(1/2*a) - 18*b*n*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) +
1/2*pi*b)*tan(3/2*a)^2*tan(1/2*a) - 18*b*n*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)
*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(1/2*a)^2 - 18*b*n*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3
/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(1/2*a)^2 + 18*b*n*x^3*e^(1/2*pi*b*
n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(1/2*a)^2
+ 18*b*n*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*
log(abs(c)))*tan(1/2*a)^2 - 18*b*n*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2
*a)*tan(1/2*a)^2 - 18*b*n*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*a)*tan(
1/2*a)^2 - b^2*n^2*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b) - 27*b^2*n^2*x^3*e^(1/2
*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b) - 27*b^2*n^2*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n
- 1/2*pi*b*sgn(c) + 1/2*pi*b) - b^2*n^2*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b) +
9*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs
(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2 + 27*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*s
gn(c) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^
2 + 27*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*lo
g(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2 + 9*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*p
i*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(
c)))^2 + 36*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*
b*log(abs(c)))*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a) + 36*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*
pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(1/2*b*n*log(abs(x)) + 1/
2*b*log(abs(c)))^2*tan(3/2*a) - 9*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*
b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(3/2*a)^2 - 27*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(
c) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(3/2*a)^2 - 27*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/
2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(3/2*a)^2 - 9*x^3*e^(
-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*t
an(3/2*a)^2 + 9*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(1/2*b*n*log(abs(x)) +
1/2*b*log(abs(c)))^2*tan(3/2*a)^2 + 27*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan
(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2 + 27*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*
b*sgn(c) + 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2 + 9*x^3*e^(-3/2*pi*b*n*sgn(x)
+ 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2 + 108*
x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)
))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(1/2*a) + 108*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/
2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(a
bs(c)))*tan(1/2*a) + 108*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(1/2*b*n*log(a
bs(x)) + 1/2*b*log(abs(c)))*tan(3/2*a)^2*tan(1/2*a) + 108*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sg
n(c) + 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(3/2*a)^2*tan(1/2*a) + 9*x^3*e^(3/2*pi*b*n*sg
n(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*a)^2 +
27*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs
(c)))^2*tan(1/2*a)^2 + 27*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*b*n*log
(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*a)^2 + 9*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/
2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*a)^2 - 9*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n
+ 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(1/2*a)^2 - 27*x^3*e^(1/2*pi*b
*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(1/2*a)
^2 - 27*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*l
og(abs(c)))^2*tan(1/2*a)^2 - 9*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(1/2*b*
n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(1/2*a)^2 + 36*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c)
- 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(3/2*a)*tan(1/2*a)^2 + 36*x^3*e^(-3/2*pi*b*n*sgn(
x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))*tan(3/2*a)*tan(1/2*
a)^2 + 9*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*a)^2*tan(1/2*a)^2 + 27*x^
3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*a)^2*tan(1/2*a)^2 + 27*x^3*e^(-1/2*p
i*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*a)^2*tan(1/2*a)^2 + 9*x^3*e^(-3/2*pi*b*n*sgn(x
) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*a)^2*tan(1/2*a)^2 - 18*b*n*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2
*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c))) - 18*b*n*x^3*e^(-3/2*pi*b*n
*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c))) - 18*b*n*x^3*e
^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c))) -
18*b*n*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*lo
g(abs(c))) - 18*b*n*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*a) - 18*b*n*x^
3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*a) - 18*b*n*x^3*e^(1/2*pi*b*n*sgn(x
) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(1/2*a) - 18*b*n*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*
pi*b*sgn(c) + 1/2*pi*b)*tan(1/2*a) + 9*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan
(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2 - 27*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2
*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2 - 27*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*s
gn(c) + 1/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2 + 9*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n -
3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2 - 9*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*
pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2 + 27*x^3*e^(1/2*pi*b*n*sgn
(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2 + 27*x^3*e^(-1/2
*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2 - 9*x
^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)
))^2 + 36*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(3/2*b*n*log(abs(x)) + 3/2*b*
log(abs(c)))*tan(3/2*a) + 36*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*b*n*
log(abs(x)) + 3/2*b*log(abs(c)))*tan(3/2*a) + 9*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*
pi*b)*tan(3/2*a)^2 - 27*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(3/2*a)^2 - 27*
x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(3/2*a)^2 + 9*x^3*e^(-3/2*pi*b*n*sgn(x
) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(3/2*a)^2 + 108*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*pi*b*n + 1/2*pi
*b*sgn(c) - 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(1/2*a) + 108*x^3*e^(-1/2*pi*b*n*sgn(x)
+ 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b)*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))*tan(1/2*a) - 9*x^3*e^(
3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b)*tan(1/2*a)^2 + 27*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2*
pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b)*tan(1/2*a)^2 + 27*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c
) + 1/2*pi*b)*tan(1/2*a)^2 - 9*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b)*tan(1/2*a)
^2 - 9*x^3*e^(3/2*pi*b*n*sgn(x) - 3/2*pi*b*n + 3/2*pi*b*sgn(c) - 3/2*pi*b) - 27*x^3*e^(1/2*pi*b*n*sgn(x) - 1/2
*pi*b*n + 1/2*pi*b*sgn(c) - 1/2*pi*b) - 27*x^3*e^(-1/2*pi*b*n*sgn(x) + 1/2*pi*b*n - 1/2*pi*b*sgn(c) + 1/2*pi*b
) - 9*x^3*e^(-3/2*pi*b*n*sgn(x) + 3/2*pi*b*n - 3/2*pi*b*sgn(c) + 3/2*pi*b))/(b^4*n^4*tan(3/2*b*n*log(abs(x)) +
3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2 + b^4*n^4*tan(3
/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2 + b^4*n^
4*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(1/2*a)^2 +
b^4*n^4*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2 + b^4*n^4*tan(1/2*b*n*log(ab
s(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2 + b^4*n^4*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^
2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2 + b^4*n^4*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(
3/2*a)^2 + b^4*n^4*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2 + b^4*n^4*tan(3/2*b*n*log(abs(x
)) + 3/2*b*log(abs(c)))^2*tan(1/2*a)^2 + b^4*n^4*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(1/2*a)^2 +
b^4*n^4*tan(3/2*a)^2*tan(1/2*a)^2 + 10*b^2*n^2*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log
(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2 + b^4*n^4*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)
))^2 + b^4*n^4*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2 + b^4*n^4*tan(3/2*a)^2 + 10*b^2*n^2*tan(3/2*b*n*
log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2 + b^4*n^4*tan(1
/2*a)^2 + 10*b^2*n^2*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)
))^2*tan(1/2*a)^2 + 10*b^2*n^2*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2 + 10*b
^2*n^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2 + b^4*n^4 + 10*b^2*n^2*tan(3/2
*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2 + 10*b^2*n^2*tan(3/2*b*
n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(3/2*a)^2 + 10*b^2*n^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*
tan(3/2*a)^2 + 10*b^2*n^2*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*a)^2 + 10*b^2*n^2*tan(1/2*b*n
*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(1/2*a)^2 + 10*b^2*n^2*tan(3/2*a)^2*tan(1/2*a)^2 + 9*tan(3/2*b*n*log(ab
s(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2 + 10*b^2
*n^2*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2 + 10*b^2*n^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^
2 + 10*b^2*n^2*tan(3/2*a)^2 + 9*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b
*log(abs(c)))^2*tan(3/2*a)^2 + 10*b^2*n^2*tan(1/2*a)^2 + 9*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(
1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(1/2*a)^2 + 9*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan
(3/2*a)^2*tan(1/2*a)^2 + 9*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(3/2*a)^2*tan(1/2*a)^2 + 10*b^2*n
^2 + 9*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2 + 9*tan(3
/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(3/2*a)^2 + 9*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2*tan(
3/2*a)^2 + 9*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2*tan(1/2*a)^2 + 9*tan(1/2*b*n*log(abs(x)) + 1/2*b*l
og(abs(c)))^2*tan(1/2*a)^2 + 9*tan(3/2*a)^2*tan(1/2*a)^2 + 9*tan(3/2*b*n*log(abs(x)) + 3/2*b*log(abs(c)))^2 +
9*tan(1/2*b*n*log(abs(x)) + 1/2*b*log(abs(c)))^2 + 9*tan(3/2*a)^2 + 9*tan(1/2*a)^2 + 9)
Mupad [B] (verification not implemented)
Time = 27.91 (sec) , antiderivative size = 122, normalized size of antiderivative = 0.76
\[
\int x^2 \cos ^3\left (a+b \log \left (c x^n\right )\right ) \, dx=\frac {x^3\,{\mathrm {e}}^{-a\,1{}\mathrm {i}}\,\frac {1}{{\left (c\,x^n\right )}^{b\,1{}\mathrm {i}}}\,3{}\mathrm {i}}{8\,b\,n+24{}\mathrm {i}}+\frac {3\,x^3\,{\mathrm {e}}^{a\,1{}\mathrm {i}}\,{\left (c\,x^n\right )}^{b\,1{}\mathrm {i}}}{24+b\,n\,8{}\mathrm {i}}+\frac {x^3\,{\mathrm {e}}^{-a\,3{}\mathrm {i}}\,\frac {1}{{\left (c\,x^n\right )}^{b\,3{}\mathrm {i}}}\,1{}\mathrm {i}}{24\,b\,n+24{}\mathrm {i}}+\frac {x^3\,{\mathrm {e}}^{a\,3{}\mathrm {i}}\,{\left (c\,x^n\right )}^{b\,3{}\mathrm {i}}}{24+b\,n\,24{}\mathrm {i}}
\]
[In]
int(x^2*cos(a + b*log(c*x^n))^3,x)
[Out]
(x^3*exp(-a*1i)/(c*x^n)^(b*1i)*3i)/(8*b*n + 24i) + (3*x^3*exp(a*1i)*(c*x^n)^(b*1i))/(b*n*8i + 24) + (x^3*exp(-
a*3i)/(c*x^n)^(b*3i)*1i)/(24*b*n + 24i) + (x^3*exp(a*3i)*(c*x^n)^(b*3i))/(b*n*24i + 24)